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LIKELIHOOD  RATIO  TESTS  FOR  RELATIONSHIPS 
BETWEEN  TWO  COVARIANCE  MATRICES 

C.  Radhakrishna  Rao 


.  1.  INTRODUCTION 

Let  S^:  p*p  and  Sg:  PXP  be  two  random  symmetric  matrices  having  Wishart 
distributions  Wp(n^,E^)  and  W^ng.Eg)  respectively,  where  n^  and  Ug  are 
degrees  of  freedom,  and  Eg  and  Eg  are  population  covariance  matrices.  In 
this  paper  the  likelihood  ratio  tests  are  derived  for  the  following  hypotheses 
on  Eg  and  Eg? 

2  2 

Hgt  Eg  »  a  £1»  °  unknown, 

2  2 

H25  ^2  *  r  +  0  ^i*  r  1-8  n.n.d.  and  p(D  «k<p,  o  unknown, 

Hj!  Eg  *  T  +  Eg,  T  is  n.n.d.  and  p(T)  •  k<p, 

H4!  ~2  *  l  +  h*  P(E)  "  k<p’ 

where  p(A)  ”  the  rank  of  the  matrix  A  and  n.n.d.  stands  for  non-negative 
definiteness. 

Applications  of  the  above  tests  to  problems  of  inference  on  "familial 
correlations"  introduced  by  the  author  (see  Rao,  1945  and  the  follow  up  in 
Rao,  1953)  are  discussed. 

The  following  well  known  results  and  notations  are  used. 

(1)  If  Eg  and  Eg  are  nonsingular,  then  Sg  and  Sg  are  nonsingular 
with  probability  1 

(il)  If  Sg  is  nonsingular,  then  there  exist  matrices  P  and  T  •  (P ')”* 


such  that 


2 


Sx  -  XT'  ,  S2  *  TAT'  (1.1) 

P'SjP  -  I  ,  p's2p  -  A  (1.2) 


where  A  is  the  diagonal  matrix  with  the  roots  of 

1  P 

f S2  —  X  S^|  »  0  as  the  diagonal  elements  (see  Rao,  1973,  p.  41). 
(ill)  Let  A:  p*p  be  a  real  symmetric  matrix  and  T^:  p*k^,  1* 1 . r 


be  such  that 

E  -  p,  p(T1:...:Tr)-p,  (1.3) 

T^A *0,  T^Tj  -  0,  i*j.  (1.4) 


Then  there  exists  a  choice  of  eigen  vectors  of 

A  such  that  the  columns  of  each  depend  on  an  exclusive  sub¬ 
set  of  R^, . . . ,Rp. 

(iv)  If  X:  p*p,  then  the  matrix  derivative  with  respect  to  X-  (X^) 
of  a  scalar  function  f(X)  is  defined  by 

<l'5> 

For  particular  choices  of  f,  we  have  (Rao,  1973, p, 72), 

Mx|  ,  ,  , 

-  |x|(X  )  (1.6) 


3  tr  MX  f 

3-x~  a  -  M  (1.7) 

3  tr  X**Hl  ,  .  , 

- = — *  -  -(X~  rl  X  )  (1.8) 

3X  "  *  " 

(v)  The  log  likelihood  of  E^,  E2  given  S^,  S2  (considering  only 
the  terms  depending  on  E^ ,  E2)  multiplied  by  2  Is 


l<!l,52'8l,82) 


-*»!  log  |  |  -  tr  E”1^  -  n2  log  |  E2 1  -  tr  e"1S2 


(1.9) 


so  that 


*31^  “”nl  ?1  +El  ?l!l  ’  l3E2J  "  ~n2  ?2  +?2?2  * 


(1.10) 


[Mote  that  in  taking  the  derivatives  we  do  not  consider  E^  and  E ^  as 
symmetric  matrices.  This  does  not  matter  so  long  as  the  optimum  solutions 
for  Z i  and  Z ^  turn  out  to  be  symmetric.] 


2.  TEST  FOR  E2  -  (o2  unknown) 

Substituting  a  Z^  for  E^  in  LCE^^IS^S^)  of  (1.9),  and  taking  deriva- 

2 

tives  with  respect  of  and  o  using  the  formulae  (1.6)-(1.8),  we  have 


+5l1(?i+0‘2?2)?i1  • 0 


(2.1) 


-  -pn.  +  o~2  tr  l71S_  ”  0. 
3o  z  ~L  ~l 


(2.2) 


From  (2.1), 


(n14n2)Z1  -  Sj^  +  o^Sj  -  T(I  +  o"2A)T' 


(n1+n2)"1l"1  -  Pd  +  o^A)"1?' 


(2.3) 


where  T  and  P  are  ae  defined  In  (1.1)  and  (1.2).  Eliminating  E^  from  (2.2) 

2 

using  (2.3),  we  obtain  the  equation  for  estimating  the  unknown  a  as 


-f  V»-! 


a  4-n  *•  5  L  7 

nlTO2  i"l  -t^+o  i“l  n2*i+ni° 


(2.4) 
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where  are  the  roots  of  |S,-A  S,  I  ■  0  and  m.  ■  n  l./n-.  The 

l  p  %1  i  1  I  z 

equation  (2.4)  has  only  one  non-negative  solution  which  we  represent  by 


a  .  Then  the  estimate  of  Z ^  is 


£  . 

,1  n1+tl2 


(2.5) 


The  likelihood  ratio  test  (LRT)  for  testing  H^:  Z^  m  a is  based  on 
the  difference 


sup  L(E1,E2|si,S2)  -  sup2  Ldj^EjJS^Sj) 
T  "  ~  ~  ~  Z^9o 


h*?2 


P  1  z 

-  log  n  [  — (  -r 

l-i  “2  -2  1  n2+Ql 
o 


n.m  +n.o  n.+n- 
2  1  1  )  1  2] 


(2.6) 


2  2 

The  statistic  (2.6)  has  an  asymptotic  x  distribution  on  (p  +p-2)/2]  degrees 
of  freedom  (d.f.)  when  n^  and  n2  tend  to  infinity. 

It  may  be  recalled  that  the  LR  test  for  E2  -  E^  (Kahlrsagar,  1978,  p.404) 


is 


01+02*1  “l4**  1 

l««  ’  »-b^>  -ir1 

1  2  m“2 


(2.7) 


which  can  be  written  as  the  sum  of  (2.6)  and 


n.+n.m.  n,+n„  _ 

m » ii1.,’1  ) 1  2 ;  ^  j. 

nl°  +n2*i 


(2.8) 


2  2 

The  statistic  (2.8)  has  x  distribution  on  1  d.f.  when  E2  ■  o  E^  and  n^,n2 

2 

are  large,  and  can  be  used  to  test  the  hypothesis  o  >1. 

2 

Further,  if  a  confidence  Interval  for  o  is  needed,  we  can  use  the 
statistic 

2  . 
no  +n_m4  n.+n  o  2n. 

log  v  [(-L— *-1)  1  2  (-j 
nl°  +n2*i 


(2.9) 
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as  x2  on  1  d.f. 

2 

A  hypothesis  of  the  type  ^  °  ^1  occurs  *n  examining  whether  two 

response  vectors  differ  by  a  scaler  multiplier.  For  example,  in  the 

evaluation  of  drugs,  two  drugs  will  be  considered  equivalent  if  their 

response  vectors,  x,y,  differ  by  a  scalar  multiplier,  since  by  a  change  of 

dosage  the  effects  may  be  made  equal.  Such  a  hypothesis  specifies  that 

2 

E(x)  *  c  E(y)  and  D(x)  ■  c  D(y).  We  have  considered  only  the  hypothesis 
2 

D(x)  -  c  D(y).  Testing  of  the  hypothesis,  E(x)  -  c  E(y) ,  under  the  con¬ 
dition  D(x)  ■  D(y)  was  considered  by  Cochran  (1943)  and  Kraft,  Olkin  and 
van  Eeeden  (1972) . 


3'.  TEST  FOR  Z2  -  T  +  o^ 


Under  the  hypothesis  r  +  c  with  T  as  an  n.n.d.  matrix  of  rank 


k  <  p  and  o  unknown,  we  can  write 


E1  "  RlRi+* * *+RpRp  "  RR* 


(3.1) 


h  "  xiRiRl+- •  *+Ak?k?k + °2<?k+i?v+i+- ' ,+RpV *RAR'  (3 * 2) 

2  2 

where  A  is  a  diagonal  matrix  with  Xj>..*>A^>o  , ...,a  as  diagonal  elements. 
We  shall  maximize 

L(El*?2lSl’?2)+tr  M^-Re')  +trM2(E2-R4R’) 

where  and  M2  are  matrices  of  Lagranglan  multipliers.  The  optimising 
equations  ara 


-  nx  Z*1  +  -  0,  -  n2  I21  +  E2ls2I21  +  “  0 


(3.3) 
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r'm,R  -  0,  i-1 . k,  l  r!m  R  -  0  (3.5) 

i-k+1 

With  o'  -  R-1,  we  have  from  (3.3) 


njl  -  u'SjU  +  R'MjR  , 

-1  i  » 

(3.6) 

n2I  -  A  U  S2U  +  R  M2R  A  . 

Adding  (3.6)  and  (3.7)  and  using  (3.4),  we  have 

(3.7) 

(nj+n2)I  -  u’SjU  +  A_1u's2U. 

The  equations  (3.6),  (3.7)  and  (3.8)  => 

(3.8) 

UiSlUi  "  ni»  UiS2Ui *  n2  V  UiSlUj  “  0  * UiS2Uj ’ 

i*  j 

(3.9) 

?i?iY  ■  °-  •  °-  . 

..;Up) 

(3.10) 

The  results  (3.9)  and  (3.10)  show,  by  using  the  results 

that  the  estimates  of  X^,  and  V  are 

(1.3)  and 

(1.4), 

*i"V  ?i““l/2!i»  1-1 . k»  !“(?lcfl;- 

where  G:  (p-k)  *  (p-k)  is  any  matrix  such  that 

..jPp)G 

(3.11) 

g’(I  +  o“2R)G  -  (n^j)  I 

(3.12) 

In  (3.12),  K  is  the  diagonal  matrix  with  the  last  (p-k)  eigen 

values  of  | ~  X  S^|  *  0  as  the  diagonal  elements.  The  equation  (3.12) 

«2  2 

together  with  (3.5)— (3. 7)  provide  the  estimate  o  of  o  as  the  nonnegative 
root  of  the  equation 


7 


(p~k)n2 

V"2 


i*k+l  n^m^+n^a 


(3.13) 


The  LRT  of  the  hypothesis  H2:  I2  -  T  +  cr Is  based  on  the  difference 


sup  L(Z1,E2|S1,S2)  -  sup  L(Z1,E2|S1,S2) 

h’-2  "  H2 


p  n  m  +n  o2  nl+n2  . 

log  n  [(-2  j-  1  )  ; 

i-k+1  nl+n2  n2  -2nl 

mi  o 


(3.14) 


which  has  a  x  distribution,  asymptotically  as  and  n2  ■+  «,  on 
f(p-k) (p-k+l)-2]/2  d.f. 

Muirhead  (1978)  obtained  a  representation  of  the  conditional  asymptotic 
distribution  of  given  and  the  last  p-k  roots  are 

equal,  and  observed  that  by  neglecting  a  linkage  factor,  this  distribution 
is  the  same  as  that  of  the  roots  of  V£V1^  where  atu*  V2  ^ave  ^ishart  distri¬ 
butions  Wp_fc(n1,B)  and  Wp_k(n2-k,o2B)  respectively.  In  such  a  case  we  may 

2 

expect  a  better  approximation  to  the  x  distribution  by  considering  the 
statistic  (3.14)  with  n2  replaced  by  n2>k. 


4.  TEST  FOR  Z2  »  r  +  Zj^  (r,  n.n.d.) 


Under  the  hypothesis  E2  •  f  +  E^  with  f  as  an  n.n.d.  matrix  of  rank  k. 


we  can  write 


.1  "  ?lh+*  *  ,+?p?p 

ll  "  Al?l?  ' *  ,+XkR  ’”^l’?k+l?kfl+'  ’  '+?p?p 


where 


(4.1) 


(4.2) 


I 
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X  >. .  .>X,  >  1. 
1  k 


(4-3) 


In  order  to  compute  the  likelihood  ratio  test,  it  is  necessary  to  ob¬ 
tain  the  estimates  of  X^,...  ,X^  subject  to  the  condition  (4.3)  which  is  some¬ 
what  difficult. 

However,  we  can  approach  the  problem  of  testing  the  hypothesis 

Z ^  *  r+ E  by  breaking  up  into  two  parts.  One  is  for  testing  the  hypothesis, 
2 

^2  *  F  +  a  (i.e.,  the  last  p-k  eigen  values  of  Z^  with  respect  to  Z^ 

2 

are  equal),  and  another  for  testing  the  hypothesis,  a  *  1  given  that  the 
first  hypothesis  holds. 

2 

The  appropriate  test  statistic  for  the  hypothesis,  Z^  9  T  +  a  is 

given  in  (3.14).  If  this  hypothesis  is. not  disproved,  ve  proceed  to  test 
2 

the  hypothesis  a  *  1  by  using  the  statistic 

_  .  nin2(p~k).l/2  rl  n  . 

Z  l2(ni+n2)  1  (<J  "l)  <44) 


which  is  asymptotically  distributed  as  a  normal  deviate.  To  deduce  the  re- 

-2 

suit  (4.5),  observe  that  a  is  a  root  of  the  equation 


(p-k)n2 


! 


nl+n2  i-k+1  i»2mi+nia2 


(4.5) 


so  that  by  the  6-method 


6  a2  -  ~  2  6m. 

p-k  i 


(4.6) 


under  the  assumption  that  the  true  values  of  , . . . ,m  are  all  equal  to  <j 
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-2 

Thus  the  asymptotic  distribution  of  a  is  the  same  as  that  of  the  average 

.  ,+mp)/(p-k)  .  Then,  using  the  results  on  the  asymptotic  distribution 
of  the  functions  of  the  roots  ^see  FanS  Krishnaiah,  1982  and 

Muirhead,  1978),  the  asymptotic  distribution  of  (4.4)  is  seen  to  be  normal. 

As  observed  earlier,  it  is  difficult  to  derive  the  LR  test  of  the 

-1 

hypothesis  that  the  last  (p-k)  roots  of  are  with  common  value 

unity.  Even  if  the  exact  LR  test  is  obtained,  it  may  not  have  an  asymptotic 
2 

X  distribution.  [It  may  be  noted  that  in  a  similar  context,  of  the  prin¬ 
cipal  component  analysis,  the  statistic  for  testing  that  the  last  (p-k) 
roots  of  a  covariance  matrix  are  all  equal  to  a  given  value,  given  by 

Anderson  (1963)  and  quoted  by  Kshlrsagar  (1978,  p.  448)  is  not  an  LRT.] 

*2 

However,  substituting  a  *  1  in  (3.14)  we  obtain  the  statistic 

P  n  +n,m  °1+n2  1 

log  n  [(  Z  J-  )  -H  (A. 7) 

i-k+l  nl  ^2  m*2 


which  provides  an  overall  measure  of  the  difference  between  the  vectors 

(m^^ . m^)  and  (1,...,1).  The  statistic  (4.7)  can  be  written  as  the  sum 

of  (3.14)  and 

p  n.+n^m.  nl+n2  2n. 

log  n  [(-1  2  — >  a  n.  (4.8) 

i»k+l  n^o  +0^1 


If  the  true  values  of  m^+^, . . . .m^  are  equal,  then  the  statistic  (4.8)  is 
asymptotically  equivalent  to 


nln2 (P"k) 
2(nx+n2) 


,-2  , .  2 

(a  -1) 


(4.9) 


which  is  the  square  of  the  statistic  (4.4),  and  hence  is  asymptotically  dis- 
2 

trlbuted  as  x  on  1  d.f.  Thus  (4.8)  is  an  alternative  statistic  to  (4.4) 
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to  test  the  hypothesis  that  the  common  value  of  the  last  (p-k)  roots  is  unity. 

In  view  of  the  remark  made  at  the  end  of  Section  3  based  on  Muirhead's 

^  2 

observation,  the  statistic  (4.7)  is  asymptotically  distributed  as  x  on 

(p-k)  (p-k+1) /2  d.f.  if  the  last  (p-k)  roots  of  are  e<lua^  with  the 

common  value  unity.  But  as  observed  earlier,  it  is  more  meaningful  to  use 
the  statistics  (3.14)  and  (4.8)  [or  (4.4)]  by  breaking  the  hypothesis  into 
two  parts,  one  specifying  the  equality  of  the  roots  and  other  specifying 
the  common  value. 

If  the  common  value  specified  is  c,  instead  of  unity,  the  statistic 
(4.8)  is  changed  to 

p  n  c+n  m  ni+n2  -2  n 

log  n  K-~~  1  )  (%-)  A].  (4.10) 

i*k+l  nia  +n2mi 

Also,  there  may  be  some  theoretical  advantage  in  replacing  n^  by 
(n^-k)  in  defining  the  statistics  (3.14),  (4.8)  and  (4.10). 


5.  TEST  FOR  I2  «  T  +  ^ 

In  Section  4,  we  considered  the  hypothesis  Z^  m  T  +  Z^  where  V  is 
n.n.d.  and  p(F)*k.  If  the  n.n.d.  condition  is  not  imposed,  then  the 
hypothesis  Z^  *  r  +  Z^  *mPl*es  that  some  (p-k)  eigen  values  of  Z^  with 
respect  to  Z ^  are  equal  to  unity.  In  such  a  case  Z^  and  Z^  can  be  written 
as 


Z,  -  R-R-+. . .+R  R 
.1  _U1  ~P~p 


ll  “  Xlh.l+'  ’  *+Xk?k?k+?k+l?k+l+‘  ‘  ,+*p?p 


(5.1) 

(5.2) 


where  need  not  be  greater  than  unity. 

The  equations  for  estimating  the  unknown  parameters  under  the  hypothesis 
E2  -  r  +  are 
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'nl?l1  +  !l1?lh1  +  h"0,  ~n2  E21  +  E21s2E21  +  M2"°  (5.3) 

M2R A  +  M^R *  0 ,  R^M^-O,  l-l,...,k  (5.4) 

where  4  is  a  diagonal  matrix  with  1,...,1  as  diagonal  elements. 

The  equations  (5.3)  and  (5.4)  are  the  same  as  those  in  (3.3)— (3.5)  except 

for  the  equation  ^  R  M  R  »  0  corresponding  to  c^. 
i-k+1  ~1~Z~1 

Proceeding  as  in  Section  3,  the  LRT  for  the  hypothesis  Z2»r  +  I^  is  seen 

to  be 


inf  log 

i 


j+(p-k)  n  m  +n  n  +n 

n  [(n"=r  >  ~r 

i-j+l  n2+nl  m.  2 


(5.5) 


2 

which  is  asymptotically  distributed  as  x  on  (p-k) (p-k+1) /2  d.f.  The  statistic 
(5.5)  is  different  from  (4.7). 

6.  FAMILIAL  CORRELATIONS 


In  an  early  paper  (Rao,  1945),  the  author  introduced  the  concept  of 
"familial  correlations"  as  a  generalization  of  the  intraclass  correlation. 
They  arose  in  a  natural  way  in  defining  a  single  measure  of  correlation 
between  members  (such  as  brothers)  of  a  family  with  respect  to  a  number  of 
measurements.  Typically  we  have  a  bxp  matrix  variable 


’xll*-",xpl 

r  1 

xi 

X  - 

.  .  .  .  . 

■ 

• 

a 

xlb’'***xpb 

(6.1) 


where  the  i-th  row  vector  corresponds  to  the  measurements  of  p  character¬ 
istics  on  the  i-th  member  of  a  family.  For  Instance,  if  we  are  considering 
brothers  in  a  family,  the  row  may  correspond  to  the  parity  of  a  brother. 

In  such  a  case,  a  natural  model  for  the  means  and  variances  and  co- 
variances  of  the  variables  In  (6.1)  Is 
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E(x^)  * 


(6.2) 


D(X) 


'  A  B  ...  B  ' 

BA  ...  B 
•  .  •  •  •  • 

,  8  B  •  •  •  A  4 


(6.3) 


Krishnaiah  and  Lee  (1974)  and  Olkin  (1973)  considered  the  problem  of  testing 
the  structure  of  the  dispersion  matrix  of  X  as  given  in  (6.3). 

In  earlier  papers  of  the  author  (Rao,  1945  and  the  follow  up  in  Rao, 
1953) ,  the  mean  vectors  were  taken  to  be  the  same  (which  is  valid  when 
the  members  of  a  family  are  not  distinguishable)  and  D(X)  is  as  in  (6.3). 
Under  this  model,  familial  correlations  were  defined  as  intraclass  correla¬ 
tions  for  suitably  chosen  linear  functions  of  the  p  measurements. 

If  we  have  observations  on  X  from  N  families,  then  we  can  write  down 
the  Analysis  of  Dispersion  (MANOVA)  for  a  two  way  classification  (families  * 
parity)  in  the  usual  way. 


Table  1.  Analysis  of  dispersion  for  two  way  classification 
by  family  and  parity 


Due  to 

D.F. 

Sums  of  squares 
and  products 
(SSP) 

Mean  squares 
and  products 
(MSP) 

E(MSP) 

Families 

N-l 

(F^) 

F 

<v 

bB  +  A-B  * 

r  +  Ex 

Parity 

b-1 

<v  ■ 

P 

<■>«> 

4  + A-B  • 

4  +  E1 

Interaction  (N-l)  (b-1)  (W^)  - 

w 

<v 

A-B  - 

h 

In  Table  1,  4  represents  the  non-centrality  parameter  which  becomes  a  null 
matrix  if  y^  are  all  equal,  and 
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ft 


F  ~  W  (N-l,  r+E.),  W  ~  W  ((N-l)  (k-l) ,  E.)  (6.4) 

p  >  -1  -  P  -1 


P  ~  Wp(N-l,  «,E1) 


(6.5) 


are  all  Independently  distributed.  The  joint  distribution  of  the  familial 

correlations  can  be  obtained  from  that  of  the  roots  of  the  equation  | F  —  XW |  *0, 
derived  by  Roy  (1939) . 

Two  hypotheses  of  interest  in  such  studies  are 


“or  Si'" "IS  ° r  ?'° 


H02!  p(!}  *  p(r)  *  k* 


(6.6) 


(6.7) 


The  hypothesis  can  be  tested  by  Wilks  A  *  |w|/|p+w|,  and  the  hypothesis 
Hq2  can  be  tested  by  using  the  statistic  (5.5),  writing  F  for  S2  and  W  for 

?1‘ 

The  familial  correlations,  P]/^****’  defined  by  Rao  (1945,  1953)  are 
the  roots  of  the  equation  |B-pA|  *0.  The  maximum  root  is  the  maximum  correla¬ 
tion  between  any  two  members  of  a  family  with  respect  to  a  linear  combination 
of  the  measurements.  The  number  of  nonzero  familial  correlations  is  equal 
to  the  rank  of  B,  and  the  hypothesis  (6.7)  is  therefore  relevant  in  drawing 
inference*  on  familial  correlations. 

It  is  seen  that  if  £^,£2» ♦ • •  are  the  roots  of  the  equation  | F  -  XW |  «  0 , 
then  the  estimate  r^  of  pi  is  obtained  from  the  relationship 


(1  +  b-l  r±) 
(b-lKl-r^  "  V 


1  ■  1,2, . . . . 


(6.8) 


If  then  the  estimates  of  *  *  ‘  *r*  obtained  from  the  roots 

.«  »• 


of  | F -  X (W+P)  |  -  0  by  the  formula 


<«-n 

n(b-l)  1-r.  *1 


(6.9) 


1 


14 


which  shows  the  relevance  of  the  hypothesis  (6.6), 

In  the  statement  of  the  hypothesis  (6.7),  no  further  condition  was 

imposed  on  B.  If  B  is  not  n.n.d.,  then  some  of  the  familial  correlations 
will  be  negative. 

Let  the  matrix  variable  X:  bxp  have  a  structure  of  the  type 


wij  +  Yj  +  €ij 


(6.10) 


where  are  constants,  and  Yj  and  are  stochastic  variables  represent¬ 
ing  family  effects  (common  to  all  members  of  a  family)  and  random  effects 
respectively,  such  that 


Then 


CovCy^e^)  -0,  CovU^.e^)  ■  0,  U*k 

D(y^» ...  ,Yp)  * ®»  * •  •  •  ™  1  ■  It...  |b.  (6. 11) 


D(X)  - 


A  B 
B  A 


B 

B 


B  B  ...  A 


(6.12) 


where  A «  B  +  C  and  B  is  n.n.d.  The  rank  of  B  is  equal  to  Che  number  of 
linearly  independent  variables  among  Yj^f-.Y  (family  effects  specific  to 
the  p  measurements).  Thus ,  if  the  covariance  matrix  of  X  is  specified  to 
be  of  the  form  (6.12),  then  the  hypothesis  (6.10)  on  the  structure  of  the 
random  variable  X  with  k  linearly  independent  y^  1*  equivalent  to 

H^:  B  is  n.n.d.  and  p(B)  -  k. 

Such  a  hypothesis  can  be  tested  by  using  the  statistics  (3.14)  and  (4.5), 
with  F  ■  S«,  W  ■  S.  and  the  corresponding  changes  in  the  degrees  of  freedom. 

«  .  A  «  .  I 

There  is  some  similarity  between  the  hypotheses  considered  in  the  present 
paper  with  those  of  Fisher  (1939)  and  Anderson  (1951).  Let  y^, . . . ,yr  be 
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the  unknovn  vectors  of  family  effects  In  N  observed  families.  The  problems 
considered  by  Fisher  and  Anderson  relate  to  hypotheses  on  y, con- 
sldered  as  fixed  parameters.  In  the  present  paper t  we  consider  y^  as 
stochastic  and  test  hypotheses  concerning  the  common  covariance  matrix  of 


I  would  like  to  thank  Professor  PtR.  Krishnaiah  for  reading  the  manuscript 
and  making  useful  comments. 
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